function z = Du(p)   
% 解析解
z(:,1) = 2^(1/2) .* pi .* sin((pi .* (4 .* p(:,1) + 1)) ./ 4)...
    .* exp(pi .* (p(:,1) + p(:,2))) .* sin(pi .* p(:,2));
z(:,2) = 2^(1/2) .* pi .* sin((pi .* (4 .* p(:,2) + 1)) ./ 4)...
    .* exp(pi .* (p(:,1) + p(:,2))) .* sin(pi .* p(:,1));
end
